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Solution of Linear Equations by an Iterative Method (SADOI Only) 

Open 

3 words at 0, 1, 2 

N words "beginning at location S3 

128 words beginning at location S5 

N words beginning at drum location S^. 

Depends on condition of equations. 

Depends on initial approximation. _ 

S3: Address of first location in memory which is used 

for latest approximation. 
Sk: Address of first location on drum which is used as 

temporary storage. 
S5: Address of first location in memory which is used 

as temporary storage. As the numbers are brought off 

the drum they are stored in these locations. 
S6: N the number of equations to be solved. 
S7: Address of subroutine for computing residuals. 
S8: Location to which control is transferred when routine 

is finished. 
135 + (Rl) = lk6 
We wish to solve the system of equation 

Ax = y 

where A is symmetric and positive definite. It is assumed 
that 



N-l 

E |a..|<l 



i = 0, 1, ..., N-l. 



This implies that the largest eigenvalue of A is < 1. 
The algorithm which is used is the following: 
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where 

r. = y - A *p . , the residual vector 
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and a is a lower bound for the smallest eigenvalue and 

b is an upper bound for the largest eigenvalue of A, 

Equation (2) is actually used during the computation in 

this program. 

Upon entry into tills routine 

contains a 

2 contains b 

A contains s n 

and S3, 1S3, ..., (N-l) S3 

contain s^^, s Q ^, ..., s^ . 



s will be made smaller during the iteration if overflow 
k 

occurs • 

Usually s, will "be a power of 2, but this is noz necessary 
and it may be more convenient to use some other base. If 
this is done, then (102) should contain the reciprocal of 
the base. The final scaling factor will be found in (506). 
In the residual subroutine , the following must be computed. 



N-l 



0, 1, ..., N-l 



This should be done in a double precision fashion using 

7^ and 70 instructions. When the subroutine is entered, 

s will be A, will contain i, and 1 will contain in 
k 

the R.H. address the address of ( s k Vfc i^' 

If there is overflow in the residual subroutine, then one 

must Jump to (100). 

This program will handle as a maximum about 600 equations. 

For a more complete discussion of the mathematical method, 

see the write up of library routine L8 - 502. 
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